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ABSTRACT 

Multi-field models of inflation predict an inequality between the amplitude tnl of the 
collapsed limit of the four-point correlator of the primordial curvature perturbation 
and the amplitude /nl of the squeezed limit of its three-point correlator. While a 
convincing detection of non-Gaussianity through the squeezed limit of the three-point 
correlator would rule out all single-field models, a robust confirmation or disproval 
of the inequality between tnl and /nl would provide crucial information about the 
validity of multi-field models of inflation. In this paper, we discuss to which extent 
future measurements of the scalc-dcpcndcnce of galaxy bias can test multi-field infla- 
tionary scenarios. The strong degeneracy between the effect of a non-vanishing /nl 
and tnl on halo bias can be broken by considering multiple tracer populations of the 
same surveyed volume. If halos down to 10 13 M & /h are resolved in a survey of vol- 
ume 25(Gpc//i) 3 , then testing multi- field models of inflation at the 3-cr level would 
require, for instance, a detection of tnl at the level of tnl ~ 10 5 given a measurement 
of a local bispectrum with amplitude /nl ~ 10. However, we find that disproving 
multi-field models of inflation with measurements of the non-Gaussian bias only will 
be very challenging, unless |/nl| ^ 80 and one can achieve a halo mass resolution of 
- 10 la M Q /h. 

Key words: cosmology: theory - large scale structure of the universe - inflation 



OO 

o 



X 



1 INTRODUCTION 



Inflation (see lLvth and Riottol (|l999h for a review) has be- 
come the dominant paradigm for understanding the initial 
conditions for the large scale structure (LSS) formation and 
for Cosmic Microwave Background anisotropy (CMB). In 
the inflationary picture, primordial densities are created 
from quantum fluctuations "redshifted" out of the hori- 
zon during an early period of superluminal expansion of 
the universe, where they are "frozen". Perturbations at 
the surface of last scattering are observable as temperature 
anisotropy in the CMB. The last and most impressive confir- 
mation of the inflationary paradigm has been recently pro- 
vided by the data of the Wilkinson Microwave Anisotropy 
Probe (WMAP) mission which has marked the beginning 
of the precision era o f the CMB measurements in space 
(|Komatsu et al.l (120111 )1. 

Despite the simplicity of the inflationary paradigm, the 
mechanism by which the cosmological curvature perturba- 
tion is generated is not yet fully established. In the single- 
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field models of inflation, the observed density perturba- 
tions are induced by fluctuations of the inflaton field it- 
self. An alternative to the stan dard scenario is rep r esente d 
by the curvaton mecha nis m dEnqvist and Sloth I (l2002h . 
iLvth and Wands I (|2002f ), iMoroi and Takahashi I (120021)') 
where the final curvature perturbations are produced from 
an initial isocurvature perturbation associated to the quan- 
tum fluctuations of a light scalar field (other than the infla- 
ton), the curvaton, whose energy density is negligible dur- 
ing inflation. The curvaton isocurvature perturbations are 
transformed into adiabatic ones when the curvaton decays 
into radiation much after the end of inflation. Alternatives 
to the curvaton model are those models characterised by the 
curvature perturbati on being generated by an in h omogene- 
ity in the decay rate (|Dvali and Gruzinov I (|2004l ') , iKofmanl 

of the particles responsible for the reheating after 
inflation. Other opportunities for generating the curvatur e 
perturbation occur at the end of inflation (iLvthl (l2005h 
ILvth and Riotto I (|2006l) ) and during preheating (|Kolb et al.l 

A precise measurement of the spectral index tiq of 
the comoving curvature perturbation £ will provide a power- 
ful constraint to single-field models of inflation which predict 
the spectral index to be close to unity. However, alterna- 
tive mechanisms, like the curvaton, also predict a value of 
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the spectral index very close to unity. Thus, even a precise 
measurement of the spectral index will not allow us to effi- 
ciently distinguish among them. Furthermore, the lack of a 
gravity-wave signal in CMB anisotropics would not give us 
any information about the perturbation generation mech- 
anism, since alternative mechanisms predict an amplitude 
of gravity waves far too small to be detectable by future 
experiments aimed at observing the £?-mode of the CMB 
polarisation. 

There is, however, a third observable which will prove 
fundamental in providing information about the mechanism 
chosen by Nature to produce the structures we see today. 
It is the deviation from a Gaussian statistics, i.e., the pres- 
ence of higher-order connected correlation functions of the 
perturbations. Indeed, a possible source of non-Gaussianity 
(NG) could be primordial in origin, being specific to a partic- 
ular mechanism for the gene ration of the cosmolo gical per- 
turbations (for a review see iBartolo et ail (|2004) ). This is 
what makes a positive detection of NG so relevant: it might 
help discriminating among competing scenarios which, oth- 
erwise, would might remain indistinguishable. 

To characterise the level of NG in the comoving cur- 
vature perturbation, one usually introduces two nonlinear 
parameters, /nl and tnl- The first one is defined in terms 
of the three-point correlator, the bispectrum, of the comov- 
ing curvature perturbation in the so-called squeezed limit 

/NL =12 pS~PS ( fc i« fc 2~ fc 3). 

fel k-2 

The second one is defined in terms of the four-point corre- 
lator, the trispectrum, in the so-called collapsed limit 



(1) 



4 Pi PS PS 



(2) 



ki fc 3 fci2 

We have normalised the correlators with respect to the 
power spectrum of the curvature perturbation, 



(3) 



and used the notation kij — (ki + kj). In all single- 
field models of inflation the bispectrum is suppressed 
in the squeezed limit and is non vanishing only when 
the spectral index deviates from uni ty , /nl = 5/12(1 
nc) — 0.02 (see Acauaviva et al. |2003j ). Maldacena 



2003) ICreminelli and Zaldarriaga I (|2004T l, ICheung et al 



20081 )1. A convincing detection of NG in the squeezed limit, 
/nl 3> 1, would therefore rule out all single-field models 
(one should be aware though that, in single-field models of 
inflation, a large NG can be generated in shapes others than 
the squeezed, e.g. in the equilateral configuration). However, 
such a detection would not rule out multi-field models of in- 
flation where the NG is seeded by light fields other than 
the inflaton. How can we derive some useful informations 
about them? In this respect, the collapsed limit of the four- 
point correlator is particularly important because, together 
with the squeezed limit of the three-point correlator, it may 
lead to the so-called Suyama-Y amaguchi (SY) inequality 
Suvama and Yamaguchi] (120081 ) , see also ISugivama et al. 1 

Based on the conditions that 



201ll ). ISmith et al. 



Eucnil (|zt 

JtaD). 



1) scalar fields are responsible for generating curvature per- 
turbations and that 2) the fluctuations in the scalar fields 
at the horizon crossing are scale invariant and Gaussian, 



Suyama and Yamaguchi proved the inequality 



tnl ^ ( q/nl 



(4) 



The condition 2) amounts to assuming that the connected 
three- and four-point correlations of the light fields vanish 
and that the NG is generated at super-horizon scales. This 
is quite a restrictive assumption. However, based on the op- 
erator product expansion, which is particularly powerful in 
characterising in their full generality the squeezed limit of 
the three-point correlator and the collapsed limit of the four- 
point correlator, it was shown that the SY inequalit y holds 
also for NG light fields (|Kehagias and Riottol (|2012f )V This 
is consequence of fundamental physical principles (like posi- 
tivity of the two-point function) and its hard violation would 
require some new non-trivial physics to be involved. 

The observation of a strong violation of the inequality 
will then have profound implications for inflationary mod- 
els. It will imply either that multi-field inflation cannot be 
responsible for generating the observed fluctuations inde- 
pendently of the details of the model, or that some new 
non-trivial (g host-like) degrees of f reedom play a role dur- 
ing inflation (|Kehagias and Riottol (|2012h ). 

Testing the SY inequality with future LSS observations 
and, therefore, the validity of multi-field inflationary mod- 
els is the subject of this paper. The squeezed limit of the 
bispectrum and the collapsed limit of the trispectrum are 
particularly interesting from the observationally point of 
view because they are associated to pronounced effects of 
NG on the clustering of dark matter halos and, in partic u- 
lar, to a strongly scale-dependent bias (|Dalal et al. I (|2008t) l. 
Measurements of the galaxy power spectrum have been ex- 
ploited to set limits on primordial non-Gaussianity c ompet- 
itive w i th those inferred from CMB o bservations (ISlosar I 
(|2009l ), iDesiacaues and Seliak I (|2010h . IXia et al. I (|201lD ). 
As we have seen, a large value of /nl in the squeezed 
limit implies that the cosmological perturbations are gen- 
erated within a multi-field model of inflation where the NG 
is sourced by light fields other than the inflaton. An in- 
escapable consequence of the SY inequality is that the 
NG is also characterised by a large trispectrum in the col- 
lapsed limit. Therefore, investigations that take advantage 
of the scale-dependent effects of NG on the clustering of dark 
matter halos should in principle take into account both /nl 
and tnl- However, since the contribution from the latter is 
suppressed by 10~ 4 (tnl//nl), setting limits on /nl under 
the assumption tnl = 0, as done in the literature, should be 
a good approximation unless |tnl| 3> /nl- 

In this paper, we will essentially try to answer the fol- 
lowing question: what values of /nl and tnl have to be 
measured in order to either confirm or disprove the SY in- 
equality ? As we shall see, even though the contributions 
from /nl and tnl are degenerate in the non- Gaussian halo 
bias, combining multiple halo mass bins can greatly help 
breaking the degeneracy. As we shall demonstrate, testing 
multi-field models of inflation at the 3-a level would require, 
for a EUCLID-like survey, a detection of a four-point correla- 
tor amplitude in the collapsed limit of the order of tnl ~ 10 5 
given a measurement of a local bispectrum at the level of 
/nl ~ 10. Conversely, we will argue that disproving multi- 
field models of inflation would require a detection of |/nl | at 
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the level of 80 or larger if dark matter halos can be resolved 
down to a mass 

The paper is organised as follows. Section 2 contains a 
short summary of the impact of primordial NG on the halo 
bias at large scales. Section 3 describes the methodology 
adopted. The last Section presents the results and discusses 
their implications. In all illustrations, the cosmology is a 
flat ACDM Universe with normalisation as = 0.807, hubble 
constant ho — 0.701 and matter content f2 m = 0.279. 



2 NON-GAUSSIAN HALO BIAS 

The effect of primordial non-Gaussianity on the halo 
bias can be computed through v a rious methods such as 
high p eaks |Matarrese and Verde I (|2008l), Shandera et al. I 
(I2011I)) or multivariate bias expan sions (|McDonald I (|200ct ). 
iGiannantonio and Porciani I (|200ct )) but, to date, the peak- 
background s plit pr o vides t he most accurate estimate 
of the e ffect (ISlosar I (120091 ) I Schmidt and Kamionkowski I 



20ld). Desiacques et al 



2011a). Smith et al. 



J2012T). 



IScoccimarro et al. I (12012^ . As shown in iDesiacques et al. I 



(|2011ah , the non-Gaussian contribution to the linear bias is 



Abi(fc) 



A N \k,z) 



{N-l)\ M s (k,z) 



(5) 



3 + 



dlnT ( a N \k, 
din o s 



where &at are Lagrangian bias parameters, S c ~ 1.68 is the 
critical threshold for (spherical) collapse and <r s is the rms 
variance of the density field at redshift z smoothed on the 
(small) scale R a of a halo. The linear matter density con- 
trast S^(z) is related to the curvature perturbation $g dur- 
ing matter domination via the Poisson equation. The latter 
can be expressed as S^(z) = M(k, z) <E>g, where 



(6) 



Here, T(k) is the matter transfer function, Q. m and Ho are 
the matter density in critical units and the Hubble rate 
today, and D(z) is the linear growth rate. M s denotes 
M(k, z)W(kRa), where W(kR 3 ) is a spherically symmet- 
ric window function (we adopt a top-hat filter throughout 
this paper). Furthermore, 



rt N Xk,z) = 



4a?P (fc) 



n 



(2tt)3 
fcjv-2, q, fc) 



M s (ki,z) 



M a {q,z) 
(7) 



is a projection factor whose fc-dependence is dictated by the 
exact shape of the Appoint function ^ JV ' of the gravita- 
tional potential. For the local constant-/NL model, the fac- 
tor J-'s is equal to /nl in the low k- limit (squeezed limit), 
so that the logarithmic derivative of J-i^' w.r.t. the rms 
variance a s of the small-scale density field vanishes on large 
scales. For all other models of primordial non-Gaussianity 
however, this term is significant for most relevant peak 
heights and becom e s negli gible in the high peak limit only 
(jDesiacaues et al71 (|2011bh ). Since the halo mass function 
may not be universal, the non-Gaussian bias correction 
should in principle be computed by taking derivative of the 



Table 1. Average host halo mass, number density, (Lagrangian) 
linear and quadratic bias factors for the low- and high-mass halo 
samples used in Fig[T] 





M(M @ /h) 


n (h 3 Mpc~ 3 ) 


6i 


62 


Halo 1 


10 12 


7 x 1CT 4 


0.2 


-0.2 


Halo 2 


10 14 


3 x 1CT 6 


2.5 


4.5 



Gauss ian halo mass function w.r.t. mass (jScoccimarro et al. I 
(|2012l )). However, because it is difficult to estimate such a 
mass derivative from real data, we will use Eq.©, which is 
valid for a universal mass function. Nevertheless, one should 
bear in mind that non-univers ality can induce additiona l 
corrections at the ~ 10% level |Scoccimarro et al. I (|2012l ); 
lMatsubara"l(|2012r i). Note also that p ath integral extensions 
of the excursion set formalism (see iMaggiore and Riotto I 
(2010)) suggest that memory terms (involving A-point cor- 
relators of the density field smoothed on any scale be- 
tween R s and Ri) co u ld also contribute a t some level 
jD'Aloisio et al. I (|2012[ ); lAshead et al. I (|2012T l). 

Specialising the above result to the bispectrum and 
trispectrum shapes considered here, the non-Gaussian bias 
correction reads 



Afei(fc) = 2/nl 



S c bi 
Ms(k) 



25 
27 



TNL 



(8) 



Ms(k) 



b 2 S c + 61 1 



dlnS 



dm a s 



where 61 , 62 are the first- and second-order Lagrangian bias 
parameters, (Jnl is another NG coefficient parametrising the 
NG arising from a cubic third-order term in the curvature 
perturbation £ and Sf^ oc (M) is the skewness of the den- 
sity field in a local, quadratic non-Gaussian model with 
/nl = 1- Strictly speaking, this expression is valid in the 
limit i C 1 only since we have ignored the fc-dependence 
of J-i 4 \ However, deviations become significant only for 
fc > 0.1 where the non-G aussian signal is negligi ble and the 
signal-to- noise saturates (|Sefusatti et al. I (|201ll )). In linear 
theory, the product o"s5^ <loc (Af) is independent of redshift. 
Therefore, at fixed values of bj and 62, the non-Gaussian 
correction induced by /nl scales as D{z)~ 1 , whereas that 
induced by <;nl and tnl does not have any extra dependence 
on redshift. For the cosmology considered here, the empirical 
relation tr s S^ 3 ' oc ~ 3.08 x 10~ 4 cr°' 145 acc urately reproduces 
the ma ss dependence of the skewness (|Desiacques et alTI 
(|2011al )). The relative amplitude of the TNL-induced scale- 
dependent bias thus is (tnl//nl)10 -4 . Whereas it is negli- 
gible in single-field inflation, it can be significant for mod- 
els with |tnl| S> ,/nl- Note that current li mits from the 
CMB trispectrum are -0.6 < t N l/10 4 < 3.3 |Smidt et al. I 
(2013)). 

An important feature of the NG bias correction is that 
its scale-dependence is degenerate in /nl, tnl and jnl in 
the large scale limit, since all the fc-dependence is then lo- 
cated in A1 s (fc)~l/fc 2 . This degeneracy can be partly broken 
by considering galaxy populations tracing halos of different 
mass and, possibly, at different redshifts. While recent stud- 
ies have analysed the problem of detecting NG through fu- 
ture large-scale surveys combining a number of observational 
datasets with simple models where only /nl is nonzero and 
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Figure 1. Confidence ellipses obtained by combining the low- 
and high- mass sample, assuming /nl = 10 and tnl = 2 x 10 4 . 



the other two nonlinear parameters are set to zero, we will 
assume here that both /nl and tnl are non-vanishing since 
we aim at testing the SY inequality ffl. We will however 
set an l to zero0. We refer the reader to lRoth and Porciani I 
(|2012h for a recent study in which both /nl and jnl are 
nonzero. Finally, it is worth mention ing that our Eq. dHI) is 
differe nt from the expression given in lGong and Yokovamal 
(2011). who neglected the mass-dependence of the skewness. 



3 METHOD 

In order to assess the ability of forthcoming experiments 
to test the SY inequality through the measurement of the 
large scale bias, we use of the Fisher information content 
on /nl and tnl from the two-point statistics of halos and 
dark matter in Fourier space. The Fisher matrix formal- 
ism has been extensively applied to predict how well galaxy 
surveys will constr ai n the nonlinear param e ter /nl (e.g., 
iDalal et al.1 (|200ct ). ICarbone etHI (|200cf ), ICunha et atl 
In particular, combining differently biased tracers 
of the same surveyed volume and weighting halos by mass 
can help mitigate the effect of cosmic variance a nd shot 
noise a nd, ther e fore, r e duce the uncertainty on /nt, jSeliak 



20111 )') 



loise a nd, thcrctorc, reduce the uncert ainty on Jnl ([s cliak 
2009 f ). ISlosar1 |2009l ). ISeliak et al. I (|2009t ). lHamaus et al. 



3.1 Fisher matrix formalism 

Here and hen c eforth , we closely follow the notation of 
lHamaus et al. I (|201ll ) and define the halo overdensity in 
Fourier space as a vector, every element corresponding to 
halos with different mass bins 



8 h = {S h (Mt),5 h (M 2 ),--- ,S h (M n )Y 



(9) 



1 Notice that this assumption also gets rid of potentia l ly lar ge 
one-loop corrections to the SY inequality l lTasinato et al.l j2012l )'). 
These corrections would be anyway below the errors we will esti- 
mate on tnl even for <?nl as large as 10 6 . 



Assuming the halos to be locally biased and stochastic trac- 
ers of the dark matter density field S, we can write the over- 
density of halos as 



S h = hS + e 



(10) 



where b is a vector whose i-component is the (Eulerian) bias 
of the i-th sample, 



bf(k,Mi,z) = l + b 1 {M i ,z) + Ab 1 (k,M i ,z) 



(11) 



and e is a residual noise-field with zero mean. We assume 
that it is uncorrelated with the dark matter. 

Computing the Fisher information requires knowledge 
of the covariance matrix of the halo samples, 



C h = (S h S n > = bb 1 P + E 



(12) 



The brackets indicate the average within a fc-shell in Fourier 
space. P = (8 2 ) is the non- linear dark matter power spec- 
trum which, on large scales, can be assumed independent of 
/nl and tnl and E = (ee T ) is the shot-noise matr i x. We 
will follow the general treatment of lHamaus et al. I (|201ll ) 
and assume that E is not simply diagonal with entries con- 
sistent with Poisson noise (see £|3.2l for explicit expressions). 

In order to simultaneously constrain /nl and tnl, it is 
pretty clear that at least two different halo samples are re- 
quired to break some of the parameter degeneracies, since 
the bias coefficients 6i, &2, the rms variance a s and the 
skewness fi£ 3 lo have distinct mas s dependences (as is ap - 
p arent from the numer ical fits of |Pe Simone et al. I (|201ll ) 
or lEnqvist et ail (|201lh ^. More precisely, the Fisher matrix 
takes the following general form 



dkk 2 1 , 



2tt 2 2 



Ti- 



dC, 



-idC h{ 



(13) 



where 0i are the parameters whose error we wish to fore- 
cast. The integral over the momenta runs from fe m i n = 
7r/(Kmrv) 1//3 to fc max = 0.1Mpc _1 //i, where V sur v is the sur- 
veyed volume and / s k y is the fraction of the sky observed. 
For illustration purposes, we will adopt the specifications 
of an EUCLID-like experiment: Kurv/sky = 25 Gpc 3 //i 3 at 
median redshift z — 0.7. We will ignore the redshift evo- 
lution and assume that all the surveyed volume is at that 
median redshift. In principle however, it should be possible 
to extract additional information on the non-Gaussian bias 
from the redshift dependence of the survey. For a single mass 
bin, the four entries of the Fisher matrix have the same k- 
dependence at low-fc. As a consequence, the determinant is 
very close to zero and, therefore, yields large (marginalised) 
errors. In this case, it is impossible to test the SY inequality 
regardless the characteristics of the halo sample, unless one 
has some prior on one of the parameters. 

In the general case of N halo populations, the entries 
of the halo covariance matrix read (o, b = 1, • • • , N) 



C ab = bf{k, M a , Z )bf(k, M„, Z )P{k) + E ab . 



(14) 



The derivative of the halo covariance matrix with respect to 
some parameter 9 is 

dC h 



(h e h' + hh{, )P, 



(15) 



where be = dh/86. We have ignored the dependence of E 
on 8 as it is expected to be small for /nl and tnl- Following 
lHamaus et al. I (|201ll ). the inverse of the covariance matrix 
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Figure 2. Halo model predictions for the 1-cr errors as a function of minimum halo mass in the limit of N S> 1 halo mass bins with 
identical number density. Red triangles show the uncertainty on /nl for a one-parameter model with /nl = (left panel) and /nl = 10 
(right panel). Filled (green) and empty (blue) squares represent the 1-cr uncertainties on /nl and tnl for the two-parameters fiducial 
models (/nLjTNl) = (0,0) (left panel) and (/nl,tnl) = (10,2 X 10 4 ) (right panel). 



Table 2. 1-tr errors obtained with N 3> 1 halo mass bins with 
M > M m i n . Top and bottom rows show results for M m ; n = 10 13 
and 10 10 M(7)/h, respectively. 





/nl = 


/nl = 10 


/nl = 


/nl = 10 




no tnl 


no tnl 


TNL = 


tnl = 2 x 10 4 


°7nl 


2.9 


3.4 


3.6 


4.2 


°"t nl 






3.6 x 10 4 


3.4 x 10 4 


"/nl 


0.5 


0.4 


1.4 


0.8 


°"t nl 






1.5 x 10 4 


0.4 x 10 4 



can be obtained using the Sh erma n-Morrison form ula (see 
ISherman and Morrisonl (|l95fj ) and lBartlettl (|l95lf )') 

r -i_ F -i E- 1 bb T E- 1 P 

Ch E ~ 1 + bTE-ibp- (16) 

On inserting the expression (|15[) into Eq. (|13p . we can write 
down the Fisher matrix for two generic parameters 0i (i — 
1,2) as 

= ^ Tr[(bb^ + b I b T )C,; 1 (bbJ + b J -b T )C^ 1 ] . (17) 

The elements of the Fisher matrix can be easily expressed 
in terms of the following quantities 

Q = bVb? (18) 

ft = b T E- 1 b ei P 

jij =hJ i B- 1 b 0] P. 

After some algebra, we obtain 

which generalises the calculation reported in lHamaus et al.~l 
(|201ll ). Note that, in what follows, 0i = /nl and 62 = tnl. 

One should bear in mind the caveat that the present 
Fisher matrix analysis assumes Gaussian uncertainties, even 



though it is likely that the estimators /nl and tnl have 
non-Gaussian distributions. One possible way of testing this 
assumption would be to generate Monte-Carlo simulations 
of the halo samples, but this is beyond the scope of this 
paper. 

3.2 Halo model predictions 

Even though the halo model makes a number of predictions 
that are not physically sensible (such as a white noise con- 
tribution in the limit k — > of the cross halo-mass power 
spectrum), it was shown to furnish a very good fit to the 
eigen values and eigenvecto rs of the halo stochasticity ma- 
trix (|Hamaus et al. I (|201Qj )). In this model, the shot-noise 
matrix can be cast into the closed form expression E, 



E = nl - bM - Mb 



(20) 



Here, M = M/p m — b(nM 2 )/2p^, M is a vector whose 
entries are the halo masses and n = n(M) is the number 
density of halos of mass M. The Poisson expectation is re- 
covered upon setting M — 0. In the limit oi N ^> 1 halo 
mass bins with identical number density n, we can replace 
the scalar products by integrals. A straightforward calcula- 
tion shows that the coefficients a, Pi and 7^ can be rewritten 
a.s 



'lol 



Pi 



(bb e 



(Mb)) 2 - {b?)(M 2 } 
)(n-J t -{Mb)) + {b 2 ){Mbe„ 



(fho\-(Mb)) -<6 2 )<M 2 ) 
(bg i b gj )n t otP 

{b 2 )(Mb ei )(Mbe 3 ) + (M 2 ){bb ei ){bbe 3 ) _ 



+ 



+ 



ntotP 

(n- \ - (Mb)) 2 - (V){M 2 ) 
((bbe % )(Mb B] ) + (bb gj )(Mb ei )) (n t -j - (Mb)) 
{n- ]-(Mb)) 2 -(b 2 )(M 2 ) 



(21) 

(22) 
(23) 

ntot-P 
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where 

(xy) = - — / dM n(M)x{M)y(M) (24) 

n tot J A/ min 

n tot = / dMn(M) = Nn . (25) 

Here, n( M) is the halo mass f u nction , which we assume to 
be of the ISheth and Tormen I (| 19991 ) form with p = 0.3, 
q = 0.73 and a normalisation A — 0.322. This yields 
(nM 2 )/pi = 75.9Mpc 3 //i 3 at redshift z = 0.7. 



4 RESULTS AND CONCLUSIONS 

We first compute the uncertainties on /nl and tnl from two 
different tracer populations and for a shot-noise matrix con- 
sistent with Poisson noise, i.e. E = diag(l/n(Mi), l/niMz))- 
We consider a nearly unbiased sample with average mass 
M ~ 10 12 : Mq /h and a high mass sample with M = 
1O 14 M //i. Table[T] summaries the characteristics of these 
populations. For a given mass M, the second-order La- 

grangian bias parameter 62 (M) is computed from the Sheth- 

(3) 

Tormen multiplicity function, whereas the skewness S s loc is 
computed from the phenomenological relation given in 
FigHJshows the resulting 68, 95 and 99% confidence contours 
for the parameters /nl and tnl when the fiducial model as- 
sumes /nl = ±10 and tnl = 2 x 10 4 . The l-cr errors are 



°7n 



23 and cr T 



2.0 x 10 b . We have tried different 



combinations of halo populations and found that the errors 
do not change significantly. At this point, we would conclude 
that galaxy bias alone cannot yield interesting constraints 
on tnl and /nl- 

The situation changes dramatically when the surveyed 
halos are divided into N ^> 1 populations of increasing 
mass, with equal number density. In Fig. [2] symbols rep- 
resent the halo model prediction for the l-cr uncertainties 
cr/ NL and ov NL in the limit of infinitely many halo bins. The 
shot-noise matrix now takes the form Eq. (|20[) . Red triangles 
indicate cr/ NL in a one-parameter model with /nl = (left 
panel) and /nl = 10 (right panel). Filled and empty squares 
represent cr/ NL and cr TNL in a two-parameters model with 
(./nl, t nl ) = (0, 0) (left panel) and (/nl, tnl) = (10, 2 x 10 4 ) 
(right panel). Results are shown as a function of the mass of 
the smallest halos resolved in the survey. Compared to the 
previous configuration, significant gains are already achieved 
for M m i n ~ 10 13 Mq/H. While the constraint on /nl is some- 
what degraded if one allows for a non-zero tnl, the l-cr un- 
certainty on tnl is of the order of (10 3 — 10 4 ), an order of 
magnitude better than in the case of two galaxy populations. 
TableHgives the l-cr errors for M min = 10 13 and 10 10 M Q /h. 

How well can we test the SY inequality with galaxy 
bias? Fig. [3] displays, as a function of /nl, the minimum 
value of tnl for which the difference (tnl — (36/25)/ NL ) 
is greater than its l-,2- and 3-cr error which, for Gaussian- 
distributed /nl and tnl, reads 



/nl» t nl 



(26) 



+ 4 



2 72 
°7 NL /nl 



1.5x10= 
1.4xl0 5 
1.3x10 s 
l.ZxlO 5 
1.1x10 s 
= 10 s 



t - 36 p > 
25 " l 


3 0" 3B,„ 


T - 36 fz > 
"25" 


2 (T T __36- ft 


T - 36 p > 
25 " L 





-50 -40 -30 -20 -10 



Figure 3. Testing the validity of the SY inequality with mea- 
surements of the non-Gaussian bias. The various curves are halo 
model predictions for N ^> 1 halo mass bins with M m ; n = 
10 13 A/q//i. At fixed value of /nl, they indicate the minimum 
Tnl required in order to have a measurement of the SY inequal- 
ity at the 1-, 2- and 3-<r confidence level. 



where a 



and 



are the entries of the in- 



°7n 



verted Fisher matrix and /nl, tnl are the values of the 
fiducial model assumed. The various curves indicate the halo 
model prediction for N ^> 1 halo populations with a min- 
imum resolved mass M m i n = 10 13 M Q //i. For instance, if a 
non-vanishing value of /nl = 10 is measured in the future, 
then the contribution induced by the collapsed limit of the 
trispectrum must be detected with an amplitude of at least 
tnl ~ 0(1) x 10 5 in order to have a 3-cr detection of the SY 
inequality with the non-Gaussian halo bias. Of course, these 
values are only indicative since the analysis is performed 
with the restrictive assumption of Gaussian errors. 

Finally, we can also assess how well galaxy bias can 
probe the violation of the SY inequality. As stated above, the 
observation of a strong violation would have profound impli- 
cations for inflationary models as it implies either that multi- 
field inflation, independently of the details of the model, 
cannot be responsible for generating the observed fluctua- 
tions, or that some new non-trivial (ghost-like) degrees of 
freedom play a role during inflation. Measuring a violation 
essentially consists in a simultaneous detection of a non-zero 
value of /nl and a (non-zero) small enough value of tnl- 
Here, we have simply estimated the smallest |/nl| such that 
(36/25)/nl is larger than the 3-cr error on tnl- Having found 
that, for the current observationally allowed range of /nl, 
the error of tnl does not significantly change if we set in all 
runs tnl = 0, we have thus computed ov NL assuming a van- 
ishing value of tnl- A comparison of 3ov NL with (36/25) / NL 
shows that, for a minimum halo mass M min = 10 13 M o /ft, 
the SY inequality cannot be tested with the non-Gaussian 
galaxy bias solely for realistic values of /nl ■ Even if halos are 
resolved down to lO lo M //i is 3cr TNL < (36/25) /nl satisfied 
only for |/nl| larger than ~ 80. 

Summarising, a large NG in the squeezed limit im- 
plies that the cosmological perturbations are generated by 
some light scalar field other than the inflaton. The SY in- 
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equality Q inevitably imposes that a large trispectrum in 
the collapsed limit is also present. However, the contribu- 
tion of tnl to the non-Gaussian halo bias is suppressed by 
10 _4 (tnl//nl) and strongly degenerate with that induced 
by /nl. Notwithstanding this, we have shown that multi- 
tracer methods can exploit the distinct mass-dependence of 
the /nl- and TNL-induced bias corrections to reduce the l-cr 
uncertainty down to ov NL < 10 4 (and simultaneously achieve 
f/ni ~ 1 — 5) for a survey covering half of the sky up to 
z « 1. The exact values depend on the mass M m i n of the 
least massive halos observed. Our results on the capability of 
testing the SY inequality through the NG scale-dependent 
bias are summarised in Fig(3] The latter shows that testing 
the SY inequality at the level of 3-<j would require detect- 
ing tnl at the level of 10 for the minimum resolved mass 
Mmin = 10 13 Mq/H. Conversely, testing the violation of the 
SY inequality requires both a much smaller resolved mass, 
Af m i n = W 10 MQ/h and a large bispectrum, [/nl| > 80. As 
mentioned above, all these results are valid provided that 
<7nl = and that the nonlinear parameters /nl and tnl 
estimated from the data are Gaussian-distributed. Relaxing 
these assumptions will be the subject of future work. 
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